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Abstract 

We present neutron scattering results for the dynamic response by superfluid and normal-fluid 
^He and the results of a simple perturbation-theory analysis which allows us to describe all aspects 
of the observed behavior by considering only density fluctuations. We show that the three key 
features that are characteristic of the dynamic reponse of superfluid ^He, as well as their dramatic 
variations with temperature, can be attributed predominantly to the Bose statistics obeyed by 
all the ^He atoms. It appears that the presence of the Bose condensate exerts at most a minor 
influence on the dynamic response. 

PACS numbers: 67.40.-w, 05.30.Jp 
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More than 6 decades after the discovery^ of superfluidity in hquid ^He at temperatures 
below T\= 2.17 K, and the almost immediate intriguing proposal that its origin lay in 
the then obscure phenomenon of Bose-Einstein condensation^, we still do not have an all- 
encompassing theory of superfluid ^He. The existence of a Bose condensate in superfluid 
■^He has since been confirmed^i^, and the current^ estimate is that about 10% of the ''He 
atoms are condensed into the zero-momentum state at T < 1 K. In a recent attemp t^i^i^° to 
develop a microscopic theory for liquid ^He, the Bose condensate has been deemed to be of 
central importance in determining the dynamic response in the superfluid phase. However, 
key predictions of this theory are inconsisten t ^^^^^ with observed behavior. 
On the other hand, the correlated basis function (CBF) approach to the superfluid state, 
which has its origin in the works of FeynmanlS and Jackson and Feenbergi^*i^ and which has 
since been developed by various group a^^' ' ^''^^^^' ' ^^' ' '^^' ' '^^ , shows that the elementary excitations 
of superfluid ^He can be calculated using a perturbative expansion based on excited states 
generated by the density operator. Good agreement between calculated and observed exci- 
tation energies has now been establishe d^'^i^'-' . In this paper, we use the CBF-theory results 
to forward a unified description of all three components in the dynamic response (elementary 
excitations, multiphonon and recoil scattering) and the variations with temperature thereof. 
Ever since the Landau's^-2. prediction of the famous phonon-roton dispersion relation for the 
elementary excitations, a large body of inelastic neutron scattering work has been accumu- 
lated (se o^i^° for references), revealing the unique features of superfluid ^He in ever increasing 
detail. Using the MARI spectrometer at the ISIS pulsed neutron facility, we have extended 
recent measurementa^^"^^^ to cover a larger dynamic and temperature range, focussing on 
the region where the excitation curve splits up into three branches. Results for the dynamic 
response functions at saturated vapor pressure (SVP) are shown in Fig. 
At the lowest temperature (r= 1.35 K, top panels of Fig. P), the spectra show three dis- 
tinct features characteristic of the superfluid phase. First, one observes the phonon-roton 
dispersion curve with the very intense roton excitations seen as the bright yellow region in 
Fig. ^ We show (a compilation2Si2£i2i of) the energies and spectral weights of these exci- 
tations in Fig. 121 Their lifetimes (related to the reciprocal of their intrinsic energy widths 
in the neutron spectra) are determined by the number of thermally excited quasiparticles^^ 
present at that temperature, which in turn govern the frequency of collisions between the 
quasiparticles. At T= 1.35 K, the observed widths primarily reflect the instrumental reso- 
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lution function (FWHM< 0.25 meV). Note also that there is no scattering at energies below 
the phonon-roton dispersion curve. 

Second, as a reddish band above and beyond the roton region in Fig. Q (and stars in Fig. 
12)), one observes a broad (in energy) response centred near the recoil energy he^ = h'^q'^/2m, 
clearly separated form the scattering on the phonon-roton dispersion curve. This component 
of the response, from which the Bose condensate fraction can be extracted-'-, mainly reflects 
the behavior of the individual ^He atoms. 

Third, there is the possibility of a neutron simultaneously exciting multiple quasiparticles. 
This multiphonon can be observed as a band of (weak) intensity located just above the 
phonon-roton dispersion curve and merging into the recoil scattering at larger q. 
In contrast, the normal- fluid dynamic response (bottom panels of Fig. ^ does not show 
any sharp excitations, nor is there any separation between the high energy recoil scattering 
and the low energy response. In addition, line-shape analysis of these and earlier— results 
indicates that the multiphonon component is absent, or at least vanishingly small in the 
normal-fluid phase. We first discuss how the features of the spectra at the lowest tempera- 
ture can be understood using (only) perturbation theory and Bose statistics. 
The shape of the phonon-roton dispersion curve is sufficient to ensure that the fluid is a 
superfluid. Specifically, the finite energy gap (roton- minumum at qr= 1.92 A~^; gap A= 
0.74 meV= 8.6 K) means that a minimum velocity (given by the slope of the dashed line 
in the top panel of Fig. |2)) is required in order to excite a roton out of the ground state, 
with frictionless superfiow for lower velocities. As Feynman pointed ouIf^, the presence of 
the energy gap is a direct consequence of the symmetry requirements placed on the wave- 
function by the Bose statistics obeyed by the ^He atoms. This energy gap also ensures that 
there can be only a few thermally excited quasiparticles at the lowest temperatures, which 
in turn leads to the long lifetime of the quasiparticles because of the vastly reduced collision 
frequency between thent^S. This explains the sharpness of the phonon-roton excitations, 
and, by consequence, also the full set of features observed in the dynamic response function, 
as we shall now demonstrate. 

Using perturbation theory, it is straightforward to derive a shape for the spectral response 
consistent with the observed spectra. Following Feynman^^, we start with a variational esti- 
mate for the excited state |g) by letting the density operator pg, which creates a density fiuc- 
tuation of wavelength 2n/q somewhere in the liquid, act on the ground state: |g) = p^|0 >. 
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The energy of this state is hepiq) = {(l\H\q) / {q\q) = h'^q'^ /2mSq, with Sq being the static 
structure factor. 

We correct this state by including the two-quasiparticle states \lm) = p^p^\0 >. In the 
following, we use the symbols \q > and \lm > to indicate that the states have been or- 
thonormalized. The corrected excited state energies are now given by the poles of the 
corresponding Green function Gi, defined by 



Given suitable expressions for the interaction matrix elements < q\H\lm > present in the 
self-energy S, the above equation can be solved in an iterative manner, and good agreement 



The validity of the above set of equations is readily demonstrated by performing a self- 



at T= K using the observed excitation energies eots as input, yielding the corrected ex- 
citation energies ep^q) — u; = eobs{q)) [which ideally would coincide with eobs{q)]- We 
plot these corrected energies, together with their corresponding intensities, as solid lines in 
Fig. 121 Clearly, the agreement with experiment is good, even near the end of the dispersion 
where the corrections to ei?(g) are very large. The discrepancies reflect inaccuracies in the 
approximations for the matrix elements, and the neglecting of states with more than two 
excitations. 

Since our primary aim in this paper is to elucidate the basic mechanisms governing the 
dynamic response of ^He instead of discussing the details of the interaction, we limit our 
discussion in the following to a simplified model, which nevertheless still contains the es- 
sential physics. Since the model reproduces the correct general shape of S, it will prove 
to be sufficient for interpreting the experimental results. In this model, we only consider 
two-quasiparticle states with momenta / and m near the roton minimum and we replace the 
matrix elements in this region by an effective Me//(g). As before, we replace the energy 
dependent corrected energies entering the self-energy S by the observed excitation energies 
heobs and lifetimes HT, yielding 



Gi = [uj- epiq) - S(g,t^)] 



-1 



(1) 



with the observed phonon-roton dispersion curve has been obtained for q < 3A U2i2ii. 



consistency check. We calculate the < q\H\lm > using a combinatio n^^i^" of convolution 
(for q < 2A~^) and superposition approximations (g > 2A^^). The self-energy is calculated 
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Determining the poles of Gi can now be done by finding the roots of the much simplified 
complex equation u + ir] = epiq) + ^{q,^^ + if]). 

We show the solution in Fig. |21 for q= 2.3 and 2.7A~^ {Mf.ff{q) was fixed by imposing 
self-consistency on the corrected excitation energies). One observes that the equation 
yields two solutions due to the appearance of a strong resonance in the selfenergy at 2 A, 
i.e., at twice the roton energy. This pronounced resonance is a direct consequence of the 
high density of states near the roton minimum in combination with the long lifetimes of 
the elementary excitations. For g-values beyond the roton (g > 2.5A^^). Eq. ^yields 
one solution with excitation energy near hepiq) (— htr in this region) and one solution 
which asymptotically approaches 2 A, almost independent of the actual values of epiq) and 
Mf.ff{q). This is exactly what is observed^^i^i in experiment (see Fig. |2I). Note that the 
linewidth of the solution near 2A is roughly given by twice the roton linewidth, as follows 
directly from Eq. |21 since the imaginary part of the denominator is given by twice the roton 
linewidth. This linewidth is still very small at T < 1 K. 

For lower g-values a similar process takes place, resulting in the lowering of the excitation 
energies with respect to tpiq) (see Fig. ^ and in the hybridization of the excitations 
making up the recoil component and the phonon-roton excitations near q= 2.3 A^^. Of 
course, the two-roton resonance disappears once q approaches 2gr (see Eq. 12)). To illustrate 
how realistic this simplified model still is, we show the predicted intensities Zq as a dashed 
curve in the lower part of in Fig. |2l The model thus shows, for all the excitations that 
constitute the phonon-roton dispersion curve, that the excitations are density fluctuations 
in character. They are not the result of knocking an atom out of the Bose condensate^*^^. 
The corrected excited state vectors are also determined in the standard mannei«2i. We 
illustrate in Fig. |3] how this gives rise to the appearance of a multiphonon component. 
The corrected single quasiparticle excited states iV'g > are obtained by rotating |g > in 
the Hilbert space spanned by [\q >, {\lm >}]. Whereas the uncorrected states \lm > were 
orthogonal to \q >, the corrected states \Xfi > are orthogonal to {ipg > (thus, < Xfi\q > ^ 
0). Since the neutron scattering cross-section of a state |.. > is proportional to | < > p, 
the multi-quasiparticle states can so be directly populated by neutrons, leading to the 
observed multiphonon component. 

We now briefly summarize the observed changes in the neutron spectra on increasing the 
temperature. There is an increasingly rapid decrease in the lifetime of the excitations^^ as 
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one approaches Tx (see Fig. [T)). This can be attributed solely to an increase in the collision 
frequency resulting from the greater number of thermally excited quasiparticles^, there 
is no need to invoke depletion of the Bose condensate. The intensity in the multiphonon 
component is transferred to the broadened single excitation^^. The separation between the 
phonon-roton dispersion curve and the recoil component first weakens and then disappears 
(see Fig. Q). The excitations near the two-roton threshold rapidly broaden (with at least 
twice the roton linewidtb^^) and weaken in intensity. 

We can understand the above changes, at least qualitatively, from the results of the simple 
model. When the lifetime of the roton decreases, the two-roton resonance will become less 
pronounced. As discussed above, it follows directly from our model that the linewidths 
of the excitations near the end of the phonon-roton dispersion curve (g > 2.5 A~^) 
increase twice as fast as the roton linewidth increases. In addition, by repeating the model 
calculations using increased roton linewidths, it follows that their intensities diminish with 
increasing linewidth so that these excitations can no longer be separated from the recoil 
componenl)^^ (see Fig. CJ. Also, with the weakening of the two-roton resonance, the excited 
state energies will be closer to hepiq)- As a consequence, the excited states \ipq > will more 
closely resemble a simple density fluctuation \q >, meaning that the multiphonon states 
|A/i > will have a smaller weight (or even no weight at all) in the neutron scattering spectra 
(see Fig. EJ. The weakening of the two-roton resonance is probably also the reason why 
the hybridization weakens (since it is the origin of the hybridization in the first place), but 
the simple model we have used here is not sufficiently accurate to allow one to follow these 
effects into the region close to the superfluid transition where the damping rates become 
comparable to the excitation energies^^. 

In summary, the results of our inelastic neutron scattering experiments and our calculations 
using perturbation theory are in good qualitative agreement. All three key features of the 
observed spectra in the superfluid phase are fully explained by the long lifetimes of the 
elementary excitations, and the changes on increasing the temperature can be attributed 
to the weakening of the two-roton resonance, the latter having its origin in the decreasing 
lifetime of the roton excitations. Our conclusion is that the dynamic response of superfluid 
''He is determined (primarily) by density fluctuations and that, although a Bose condensate 
exists in superfluid ^He, it plays at most a minor role in determining this response. 
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FIG. 1: The (resolution broadened) dynamic response of liquid ^He at SVP for five temperatures 
(T=1.35, 1.85, 2.0, 2.1 and 2.3 K, from top to bottom). We also obtained results (not shown) at 
2.15, Tx-O.OOS and 3.0 K. The panels on the right show expanded views of the roton region. 
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FIG. 2: Top: the phonon-roton dispersion curve of superfluid ^He at T < 1.3 K (circles). The 
roton excitation is given by the minimum of the dispersion curve at qr= 1.92 A~^, with energy gap 
A. For g-values beyond the roton, the dispersion curve asymptotically approaches 2A (dotted line) 
and terminates at about 2qr. The stars indicate the peak positions of the recoil component of the 
scattering. The uncorrected excited state energies hep{q) are given by the dashed-dotted curve, and 
the corrected energies (self-consistent calculation) by the solid line. Bottom: corresponding weights 
Zq of the phonon-roton excitations (circles: experiment, solid line: self-consistent calculation, 
dashed-dotted line: simple model). 
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FIG. 3: Solution of the roots oi ui + ir] = epiq) + S(g,a; + irj) for q= 2.3 (solid line) and 2.7 A'^ 
(dashed line). Only the real part of the equation is depicted here. In each case there are two 
solutions, shown by circles. The apparent third solution does not satisfy the imaginary part of the 
equation. 



11 




FIG. 4: Pictorial illustration of how an excited (single quasiparticle) state \ipq >, which does 

not coincide witli a simple density fluctuation \q >, necessarily leads to the observation of multi- 
quasiparticle states {X/j, >. 
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